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(N : ABSTRACT 

^ . 

We present a Very Long Baseline Interferometry image of the water maser 
^ ■ emission in the nuclear region of NGC 3393. The maser emission has a linear 

distribution oriented at a position angle of ~ —34°, perpendicular to both the 
i-^' kpc-scale radio jet and the axis of the narrow line region. The position-velocity 

I ■ diagram displays a red-blue asymmetry about the systemic velocity and the esti- 

^ • mated dynamical center, and is thus consistent with rotation. Assuming Keple- 

^ ■ rian rotation in an edge-on disk, we obtain an enclosed mass of (3.1±0.2) x 10'' Mq 

within 0.36 ± 0.02 pc (1.48 ± 0.06 mas), which corresponds to a mean mass den- 
sity of ~ IO^'^MqPc"^. We also report the measurement with the Green Bank 
Telescope of a velocity drift, a manifestation of centripetal acceleration within 

■ the disk, of 5 ± 1 km s ~^ yr~^ in the ~ 3880 km s~^ maser feature, which is 

■ most likely located along the line of sight to the dynamical center of the system. 
From the acceleration of this feature, we estimate a disk radius of 0.17 ± 0.02 
pc, which is smaller than the inner disk radius (0.36 ± 0.02 pc) of emission that 

O ■ occurs along the midline (i.e., the line of nodes). The emission along the line of 

sight to the dynamical center evidently occurs much closer to the center than the 
^ . emission from the disk midline, contrary to the situation in the archetypal maser 

systems NGC 4258 and NGC1068. The outer radius of the disk as traced by the 
masers along the midline is about 1.5 pc. 

Subject headings: galaxies: active — galaxies: individual (NGC 3393) — galaxies: 
Seyfert — ISM: molecules — ISM: jets and outflows — masers 
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1. Introduction 



^NGC 3393 is a face-on early spiral (Sa) galaxy classified optically as a Seyfert 2 (IVeron-Cetty &: Veron 

20031 ). It displays a prominent ring and interacts weakly with a companion 60kpc away 
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( ISchmitt et al.ll2001al ). There is compelling evidence that the galaxy contains an active galac- 
tic nucleus (AGN). Its X-r ay spectrum exhibits a 6.4 keV Fe Ka line ^ is an unambiguous in - 
dicator of nucle a r acti vity (IMaiolino et al.lll998l : lGuainazzi et al.ll2005l : iLevenson et al.ll2006l ). 



Levenson et al.l (120061 ) estimate the AGN's 2 — 10 keV intrinsic luminosity to be 10^'^ Lq based 



on modelling of the Fe line in the Chandra X-ray Observatory data. A Compton-thick, cold 
reflection model applied to the BeppoSax X-ray data yields an obser ved 2 — lOkeV flux of 
3.9 X 10~^^ergcm^^ s^^ absorbed by a c olumn density of > 10 



25cm-2 (IMaiolino et al.lll998[ ). 



If we assume, following IComastril (120041 ). that the observed 2 — lOkeV flux is 1 — 10% of the 
unabsorbed 2 — lOkeV flux due to reflection and scattering, then the intrinsic 2 — lOkeV 
luminosity is IO^'^'^'^Lq, which is consistent with the luminos i ty det ermined from the Fe 
Ka line. Based on XMM and BeppoSax data, iGuainazzi et al. J2005h report an absorbing 
column density of (4.4;^^;^) x 10^^ cm^^, a spectral or photon of index 2.8~^qj, and an observed 
flux of {QIL^) X 10~^'^ergcm~^ s~^, which corresponds to an intrinsic 2 — lOkeV luminosity 
Lq (as computed by the Portable Interactive Multi-Mission Simulatoi0). In the 
context of the AGN unifled model, such a high column density is indicative of an almost 



edge-on obscuri ng structure along the line of sight to the nucleus ( iLawrence fc Elvislll982 



Antonuccil 119931 ) . Furthermore, images of the inner- kpc region reveal outflows from the nu- 
cleus, consistent with the presence of an AGN. The Narrow Line Region (NLR) as traced by 
[O III] emission has an S-shaped morphology with an opening angle of ^ 90° and extends 

o i 

^ 40 PC on either side o f the nucleus along position angle (P.A.) of ~ 55 (ISchmitt &: Kinney 



19961 : ICooke et al.ll2000l ). Very Large Array (VLA) and Australia Telescope Compact Ar- 
ray (ATCA ) observations reveal a double- sided jet with a total exte nt of ~ 700 pc along 
P.A.~ 56° JSchmitt et al.lboOlbh or P.A.~ 45° jMorganti et al.lll999[ ). The existence of an 



AGN in the nucleus of NGC 3393 is thus well established based on the available data. 

NGC3393 does not show evidence of signiflcant nuclear star formation. The IR AS 
fluxes measured on an ~ 30 kpc scale (assuming IRAS beam of ~ 2'; iMoshir et al.lll990l ) are 



consistent with a total infrared luminosity of 10 Lf^, a dust mass o f ~ 5 x 10 Mq, and two 



relatively cool dust components at 30 and 130 K (jPiaz et al.lll988l ). Under the assumption 
that star formation alone is responsible for the infrare d flux, the IRAS lu minosity yields 
a star formation rate of ~ 4 MQyr^^ on kpc scales (IVeilleux et al.l 1 19941 ). which is not 
consiste nt with a high level of star -forming activity. Detailed stellar population synthesis 
work by ICid Fernandes et al. J2004h conflrms that the blue optical spectrum of the central 
~ 200 pc is dominated by an old stellar population (> 2.5 x 10^ yrs) with a small contribution 
(14%) to the featureless continuum from an AGN. Hence, starburst activity probably does 
not play a signiflcant role in the nucleus of NGC 3393. 



"'^ http://heasarc.gsfc. nasa.gov/Tools/w3pinims. html 
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The NGC 3393 nucleus is also a source of water maser emission, which is currently the 
only resolvable tracer of warm dense molecular gas in the inner parsec of any AGN beyond 
~ 1 Mpc. The maser spectrum of NGC 3393 shows a characteristic spectral signature of 
rotation in an edge-on disk: two complexes (~ 70mJy) symmetrically offset by ~ 600kms~^ 
from the systemic velocity (henceforth, high-velocity emission) and a single spectral complex 
(~ 28mJy) within 13 0kms~^ of the systemic velocity (henceforth, low- velocity emission; 



Kondratko et al.ll2006l ). Very Long Baseline Interferometry (VLBI) maps of seven water 



maser sources that show similar spectral signatures have been interpreted in a context of a 



to 1 pc from a supermassive black h( 
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The maser spectrum of NGC 3393 thus provides indirect evidence for an edge-on pc-scale 
molecular disk orbiting a supermassive black hole. 

VLBI maps of nuclear water maser emission have been used to accurately estimate 
black hole masses and pc-scale accret ion disk geornetries . In three w ater maser systems 



mapp ed with the VLBI — NGC 4258 flMivoshi et al 



19971 ). and the Circinus Galaxy (iGreenhill et al.ll2003l ) — position and line-of-sight velocity 



19951 ). NGC 1068 f lGreenhill fc Gwinn 



data provided evidence for differential rotation and enabled accurate estimation of black hole 
mass and pc-scale molecular disk structure. In another system, NGC 3079, the rotation curve 
traced by the maser emission appears flat and was interpreted in the cont ext of a pc-scale. 



thick, edge-on, self-gravitating, and possibly star forming molecular disk (IKondratko et al. 



20051 ). In addition to mapping pc-scale molecular disk structure and accurately determin- 
ing the mass of supermassive black holes, nuclear water maser emission has also been used 
as a distance indicator. Distance determination is possible for systems where a detailed 
knowledge of the structure of the disk from VLBI is combined with a measurement of either 
maser proper motions or drifts in line-of-sight velocity of spectral features (i.e., centripetal 
acceleration). The distance to NGC 4258 obtained in this manner is the most accurate ex- 
tragalactic distance thu s far, is independent of standard candle calibrators such as Cepheids 
( iHerrnstein et al. I ll999h . and ha s contributed to the analysis of the Ceph eid period-luminosity 
relation (Freedman et al. 2001: iNewman et al.ll200ll : iMacri et al.ll2007l ). 



In this work, we present a VLBI map of the pc-scale accretion disk as traced by water 
maser emission and estimate the mass of the black hole in NGC 3393. Data calibration 
and reduction techniques are discussed in Section [2J Spectral-line images of the inner-pc 
region are presented in Section [31 In Section |H we interpret the observed kinematics of the 
maser distribution in terms of a nearly edge-on pc-scale molecular accretion disk that orbits 
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a central mass. In this work, we a dopt systemic velocity for NGC 3393 based on 21 cm line 
measurements of 3750 ± 5kms~^ ( iTheureau et al.lll998l ). which yields a distance of 50Mpc 
for a Hubble constant of 75kms~^ Mpc~^. With this adopted distance, Imas corresponds 
to 0.24pc. 



2. Observations and Calibration 

NGC 3393 was observed in spectral-line mode with the Very Long Baseline Array (VLBA) 
of the NRAol, augmented by the phased VLA and the Green Bank Telescope (GBT), for 
approximately 5 hours on 2004 April 26/27. The source was observed with 8 x 16 MHz 
intermediate-frequency (IF) bands that covered an uninterrupted velocity range of 3048 to 
4430 km s~^. Limits on the VLBA correlator data rate necessitated two processing passes 
through the correlator (with 256 channels per IF) and resulted in two separate data sets 
that do not overlap in velocity: one in left-circular polarization covering the velocity range 
from 3048 to 3737 km s~^, the other in right-circular polarization extending from 3737 to 
4430 km s~^. The data were reduced using the Astronomical Image Processing System. The 
amplitude calibration included corrections for atmospheric opacity. Antenna gain curves and 
measurements of system temperature were used to calibrate the amplitude data for each of 
the VLBA stations. Amplitude calibration for the VLA was based on the measurements with 
the VLA of flux densities for VLBI calibrators with respect to 3C286, for which we adopted 
a 22-GHz flux density of 2.56 Jy. To minimize systematic errors due to atmospheric opacity 
effects, we used in this calibration only VLA scans of VLBI calibrators with elevations close 
to the elevation of 3C286 {elev = 28°, 29°, 54°, 43°, 64° for 3C286, 1055-248, 4C 39.25, 
0727-115, 0748+126, respectively). 

The data in each polarization were corrected for the effect of the parallactic angle of 
the source on fringe phase. The zenith path delay over each antenna was estimated via ob- 
servations of 10 strong (> 0.5 Jy) calibrators with sub-mas positions for ~ Ihour before and 
~ Ihour after the main imaging experiment. The residual multiband delays after applying 
this calibration were < 0.2 ns, which corresponds to a phase error of less than 1° for a 16 MHz 
IF band (see Tabled]). Time variation in delays and fringe rates due to the troposphere and 
clock uncertainties as well as electronic phase difference among bands were removed via 
~ 4min observations of strong (3 — 7Jy) calibrators (0727-115, 0748+126, 4C 39.25) made 
approximately every hour. We estimate that this calibration is accurate to within 0.2 ns and 



^The National Radio Astronomy Observatory is operated by Associated Universities, Inc., under cooper- 
ative agreement with the National Science Foundation 
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1 mHz for delays and rates, respectively; we note that a residual rate of 1 mHz results in a 
phase error of 0.5° (if not corrected) over a 1 hour time span, the approximate time separa- 
tion between the calibrator scans (Tabled]). The complex bandpass shapes of each antenna 
were corrected using the data on 0727-115, 0748+126, and 4C 39.25. Phase and amplitude 
fluctuations due to the troposphere and clock uncertainties were removed via observations 
every ~ 100s of 1055-248, located 2.4° away from NGC3393, self-calibrating the data on 
1055-248, and applying the resulting phase and amplitude solutions to the target source. 
As a result of this process, the position of the maser on the sky was determined relative 
to 1055-248. Based on VLA data obtained at two epochs in B and BnA configuration, we 
estimate the 22 GHz flux of 1055-248 to be 0.38 ± 0.03 Jy. This measurement as well as the 
photometric data from the NASA extragalactic database (NED) and from NRAO VLA Sky 
Survey (0.33 Jy at 1.4 GHz; ICondon et al.lll998l ) suggests that 1055-248 is a fiat spectrum 
source and is therefore expected to be unresolved on mas scales. In fact, using the source 
model obtained from self-calibration, we estimate that 90% of 1055-248 flux originates from 
a point source. The calibrator appears nearly point-like also in a narrow bandwidth (250 kHz 
~ 3.4 km s^^), which confirms the quality of bandpass calibration. Based on the considera- 
tion of the tropospheric str ucture function as predicted by the Kolmogorov turbulence theory 
( ICarilli &: Holdawaylll999l ). we estimate that the interferometer phase towards 1055-248 re- 
produces that towards NGC3393 to within ^rms < 40°, (assuming a water vapor scale height 

o 

of ~ 2 km, a representative antenna elevation angle of ~ 25 , and a tropospheric pattern 
speed of < 15ms~^; Tabled]) which corresponds to a loss of signal on the target source of 
at most 1 — exp[— 0^jjjg/2] = 0.22. The uncertainty in phase due to calibrator-target angular 
separation is the most significant source of error in this experiment and depends sensitively 
on the assumed pattern speed, which is not well constrained for this experiment (Tabled])- 

After calibrating and imaging the spectral-line data set, the positions and the integrated 
fluxes of the maser features were obtained by fitting two-dimensional elliptical Gaussians to 
the distribution of the maser emission on the sky in each spectral channel. The resulting 
uncertainties in maser spot positions, based on formal error analysis, are consistent with the 
theoretical expectations based on signal-to-noise ratio (SNR) and beam size. The synthesized 
beam dimensions were 3.5 x 1 mas at P.A.~ —1° and the resulting noise level was ~ 5 mJy 
in the ~ 3.4kms~^ spectral channels, a spectral breadth comparable to the linewidths of 
the maser features. To estimate the fraction of the power imaged by the interferometer and 
to monitor the maser with the intent of measuring velocity drifts, we acquired single-dish 
spectra of NGC 3 393 w ith the GBT using the configuration and calibration described in 



Kondratko et al.l (120081 ). All velocities have been computed in accordance with the optical 



definition of Doppler shift and are in the heliocentric reference frame. 
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3. Results 

The spectrum of imaged power agrees to within la with the total power spectrum 
obtained with the GBT about 9 months after the VLBl observation (Fig. [1] and Table [2]). 
The difference between the two spectra may be due to source variability, since single-dish 
monitoring of water maser sources has reveale d substantial flux variabi lity on timescales of 



months to years (this work for NGC3393 and iBaan fc Haschickl (119961 ) for other sources). 
The maser emission is distributed on the sky in a linear arrangement at P.A.~ —34° (Fig. 
|2] and Table [2]), perpendicular to both the kpc-scale radio jet (P.A.~ 4 5°, Morganti et al. 



1999 ; ~ 56° Schmitt et al. 2001b) and the axis of the NLR (P.A.~ 55°; ISchmitt fc Kinnev 



19961 : ICooke et al.ll2000l ). The maser emission on the sky is clearly systematically distributed 
according to velocity: emission red- and blue-shifted with respect to the systemic velocity are 
located in the north-western and south-eastern parts of the image, respectively. Although 
the two VLBA correlator passes did not overlap in velocity, we can nevertheless place a limit 
on their registration by aligning the common low velocity maser spots in the two data sets. If 
we assume that the mapped low-velocity features arise at the same location on the sky, then 
the two correlator passes are registered to within 0.3 ±0.1 mas and 0.9 ±0.3 mas in RA and 
DEC, respectively. The systematic offset between the results of the two correlator passes was 
also estimated by applying the calibration of one data set to the other. The transfer of phase 
and amplitude solutions from self-calibration on 1055-248 from one correlator pass to the 
other resulted in an offset of 1055-248 from map center of 0.30±0.02 mas and 0.52±0.06 mas 
in RA and DEC, respectively, which is consistent with the offsets based on the location of the 
low-velocity maser features. These offsets provide an estimate of the systematic uncertainty 
in the registration of the two velocity sections of the maps. Note that if we were to apply 
these offsets, the high velocity masers would better fit a straight line distribution, and the 
low velocity features would be more tightly clustered (see Figure [2]). However, this shift has 
an insignificant effect on the disk parameters derived in the next section. Radio maps of 
the inner-kpc are suggestive of a jet perpendicular to the distribution of maser emission on 
the sky. However, we find no evidence in our data at a 4(T = 2.4 mJy level for the 22 GHz 
continuum on ~ 50 pc scale. 

The low-velocity spectral features are expected to drift in velocity since they are pre- 
sumably moving across the line of sight where the centripetal acceleration vector is along the 
line of sight, and therefore at a maximum. To determine these velocity drifts, we monitored 
NGC 3393 with the GBT and report the measurement of the centripetal acceleration in the 
low- velocity comp l ex at ~ 3880 km s~^ (Fig. [3] and Table [3]). Using the code described in 



Humphreys et al.l (120081 ). we decomposed the spectra for the six epochs into three Gaussian 
components and used an iterative least squares technique to solve simultaneously for velocity 
drifts, amplitudes, peak centroids, and widths of all Gaussian components at all epochs. We 
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obtain a reduced of 1.2 and velocity drifts of 5.3 ± 0.7, 5.3 ± 0.2, and 4.4 ± 0.2 km yr~^ 
for components at 3871.0 ± 0.9, 3874.9 ± 0.2, and 3879 ± 0.2 km s respectively, where 
the reference epoch for component velocities is 2005 January 15. In addition, spectra ob- 
tained since 2004 October 17 reveal a strong (< 0.3 Jy) feature at ~ 4051 km s~^ that was 
present neither in earlier spectra of the source nor in the VLBI map (Fig. [3] and Table [3]). 
The first three epochs on this feature showed a negative centripetal acceleration of about 
— 4kms~^ yr~^, as would be expected for low- velocity emission that arises from behind the 
dynamical center; however, the line stopped drifting in the more recent spectra. We note 
that such behavior can be reproduced by a variation in strength of multiple components that 
do not drift in velocity (i.e., stationary components). In fact, from a Gaussian decomposition 
of the five available epochs, we infer that the data on this complex are consistent (reduced 
of 2.0) with three time variable but stationary (< 0.7kms~^ yr"^) Gaussian components at 
4051, 4052, and 4053 km s'^ High-velocity emission is not expected to drift in velocity since 
it is located along the disk midline where the centripetal acceleration vector is perpendicular 
to the line of sight. We thus suggest that the newly detected complex at ~ 4051 kms~^ may 
be a high- velocity emission component. We note that low- velocity emission that arises from 
behind the dynamical center has not been detected to date from any known nuclear water 
maser sources. This non-detection can be explained if fre e- free absorption by an interven- 



ing ionized gas is considerable (e.g., iHerrnstein et al.lll996l ) or the presence of a background 



nuclear continuum source is necessary to generate, via ampli fication, low-velocity emission 



luminous enough for us to detect (e.g., IHerrnstein et al.lll997l ) 



4. Discussion 

We interpret the linear distribution of the maser emission, perpendicular to the radio 
jet and to the axis of the NLR, as well as the segregation of the blue- and the red-shifted 
emission on the sky in the context of a pc-scale molecular disk. We assume that the disk is 
nearly edge-on and, based on the measured positive centripetal acceleration, we infer that 
the low-velocity emission lies in front of and along the line of sight to the dynamical center. 
The distribution of maser emission on the sky is consistent with a relatively straight (i.e., 
nonwarped) disk. We note that a 4th degree polynomial fit to the maser distribution yields 
a marginal (i.e., ~ 10%) improvement in over a straight line fit. Thus, the evidence for a 
warp is at most tentative with our SNR. 

We use the mean position of the low-velocity maser features to estimate the location of 
the dynamical center (Fig. [2]), RA^j/ and DECsj^. The absolute position of the dynamical 
center is 
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^2000 = 10^^48^23^4659 ± 0^0001 
S.nnn = -25°09'43M77 ± O'.'OOl 



The errors are dom inated by the uncertainty in the position of 1055-248 of about 1 mas 
jBeaslev et aLll2002h . 



The impact parameter of each maser feature was calculated as 
a/ {RAi — RAbhY + {DECi — DECbhY- The resulting position-velocity diagram displays 
a red-blue antisymmetry about the adopted Vsys and estimated dynamical center and is thus 
consistent with rotation (Fig. H]). From a fit of the Keplerian rotation to the high- velocity 
features (x/j = 0.6), we obtain a mass of (3.1 ± 0.2) x 10^ Mq enclosed within 0.36 ± 0.01 pc 
(1.48±0.06 mas; to estimate the minimum impact parameter from high-velocity data, we used 
the probability distribution function for a minimum of random variables following the Rice 
distribution) . The fit to the blue- or the red-shifted emission alone yields (3.0 ±0.3) x 10^ Mq 
or (3.2 ±0.2) X 10'' M©, respectively. If one correlator pass is shifted 0.9 mas so that its single 
low-velocity maser feature overlaps the two low-velocity features in the other data set, then 
the mass becomes (3.5 ± 0.2) x W Mq enclosed within 0.41 ± 0.02 pc (1.71 ± 0.07mas), 
which reflects the sensitivity of our results to systematic errors. We note that, in addition to 
the Keplerian rotation, the data are also consistent (in the sense that xlj ^ 1) with f oc r" 
for — 1 < a < —0.1. In particular, we obtain a minimum of = 0.4 at a = —0.2, which 
suggests that the disk might have significant mass with respect to that of the black hole. 
The central mass would be 2.6 x 10^ enclosed within 0.36 pc (if we assume spherical 
symmetry). The mass of the disk traced by the high- velocity maser emission (0.36 — 1 pc) 
is 1.9 X 10^ M0, which can be compared to disk r nasses (computed on s imilar scales) of 
~ 7 X 10^ and ~ 9 x 10^ M© for NGC3079 jKondratko et all l2005h and NGC1068 



(ILodato fc Bertirul2003l ). respectively. A flat rotation curve model, which gives 1.1 x 10^ Mq 
enclosed within 0.16 pc, is excluded by the data {x\ ~ 20), unless velocity dispersion on 
the order of 30kms~^ is included. Such a large velocity dispersion would most likely be 
indicative of macroscopic random motions among the molecular clumps responsible for the 
maser emission rather than turbulence within the clumps. 

The mean mass density corresponding to 3.1 x 10^ Mq enclosed within 0.36 pc is 
~ 10^'^ Mq pc~^. The relatively high mean mass density for NGC3393 is suggestive of 
a massive central black hole, which is consistent with the X-ray observations of the nu- 
cleus. The estimated enclosed mass of (3.1 ± 0.2) x 10^ Mq is in agreement with the 
empirical relation between bulge velocity d i spersion and b l ack h ole mass (Gebhardt et 
al. 2000a, 2000b; iFerrarese k MerrittI I2OO0I : iFerrarese et aD boOlh . If we adopt Mbh = 
1.2 X 10^ Mq [(t/(200 km s"^)]3-^^ (Gebhardt et al. 2000a, 2000b), then the velocity dis- 
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persion of the bulge, 184 ± 18kms ^ from the central ~ 35kpc ( Terlevich et al.lll990l ) or 
157 ± 20kms~^ from the central ~ 200 pc (ICid Fernandes et al.ll2004l ). predicts black hole 



masses of (9 ± 4) x 10^ Mq or (5 ± 2) x 10^ Mq, respectively, roughly a factor of two higher 
than our measurement. 

The most reliable estimate of the AGN's 2 — lOkeV int rinsic luminosity is 10^'^ Lq de- 



termined from Fe Ka line l uminosity (ILevensqn et a. 



spectral energy distribution (IFabian fc Iwasawa 



1999 



20061) ■ Based on an average quasar 



Elvis et al.ll2002l ) , the 2 - 10 keV lumi- 



nosity is 1 — 3% of the AGN bolometric luminosity and we obtain a bolometric luminosity for 
NGC 3393 of lO^o-^-^^ ^ Lq. We note that this estimate for the AGN bolometric luminosity is 
consis tent with the total IRAS luminosity of the source of ~ 10^° L© (measured on ~ 30kpc 



scale; iMoshir et al.lll990l ). The Eddington luminosity of a 3 x 10" M© object is 10 L 



0- 



Assuming that all of the enclosed mass is concentrated in a supermassive black hole, the 
bolometric luminosity of the central engine yields an Eddington ratio of 0.06 — 0.2, which is 
consistent with the 0.01 — 1 range obtained for Seyfert 1 galax ies, representative supermassive 



black hole systems (e.g., |Padovanilll989l : IWandel et al. 



19991). but larger than for adyectio n 



Herrnstein et al.lll998l : lYuan et al.ll2002l ). 



dominated accretion flow systems (e.g., io~^-^^^ 
Assuming a standard accretion efficiency of ~ 0.1 (Frank, King, & Raine 2002; see also 
Marconi et al. 2004), we estimate a mass accretion rate of M = 0.04 — 0.1 Mq year~^. 

For a central mass of (3.1 ± 0.2) x 10^ M© and the measured centripetal acceleration of 
a = 5 ± 1 km s~^ yr~^ , we estimate disk radius of the systemic maser feature at ~ 3880 km s~^ 
^GMsH/a = 0.17±0.02pc, which is significantly smaller than the inner disk radius 



of r 



sys 



of the high- velocity emission (0.36 ±0.02pc). Evidently, the systemic emission in NGC 3393 
arises much closer to the dynamical center than the high-v elocity emission, which is in con 



trast to the situation in NGC 4258 (IHerrnstein et al.ll2005[ ) and NGC 1068 fIGreenhill et al. 



19971 ). where disk radii of low- velocity features is about equal to the inner radius of the high 
velocity masers. It has been suggested that the systemic er nission in NGC 4258 re sides in 
a bowl that is a consequence of an inclination- warped disk (IHerrnstein et al.ll2005l ). Such 
a warp in the accretion disk structure might also determine the preferred radial location of 
the low- velocity features in NGC 3393. The resulting orbital velocity of the ~ 3880 km s~^ 
maser feature is 890 ±60 km s~^ (which might be as high as 920 ±60 km s~^ due to systematic 
errors). We note that the 130kms~^ offset of this feature from the adopted systemic velocity 
might be due its location within the disk at a non-zero azimuthal angle from the line of 
sight to the central engine. Using this velocity offset, we estimate = 8° and the resulting 
corrections to the derived radius and orbital velocity are much smaller than the correspond- 
ing uncertainties. The newly detected feature at ~ 4051 km s~^ that we postulate to be 
high-velocity emission was not detected in the VLBI experiment but would appear at a large 
disk radius which from the computed enclosed mass and Keplerian formula, r = GMbh/v'^, 
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is ~ 6 mas or 1.5 pc (see Figure H]). Hence, we estimate that the accretion disk extends from 
0.17pc to 1.5 pc. 



5. Conclusion 

We have mapped for the first time the maser emission in the nuclear region of NGC 3393. 
We interpret the linear distribution of the maser emission and the segregation of the blue- 
and the red-shifted emission on the sky in the context of a pc-scale nearly edge-on molecular 
disk that orbits a central mass of (3.1 ±0.2) x W Mq enclosed within 0.36 ± 0.02 pc (1.48 ± 
0.06 mas). We also report the measurement of centripetal acceleration, a = 5±1 kms^^ yr^^, 
in the low- velocity maser feature at ~ 3880 km s~^, which yields disk radius of 0.17±0.02pc 
for the derived central mass. The low- velocity emission in NGC 3393 occurs much closer 
to the dynamical center than the high-velocity emission, in contrast to the situation in 
NGC 4258 and NGC 1068, two archetypal maser systems. An independent estimate for the 
disk radius of the low-velocity features would be provided by the measurement of their 
proper motions. For a distance D, a central mass 3.1 (-D/50Mpc) x 10^ Mq, and a radius 
0.17 (D/50 Mpc)^/^ pc, we expect motions of ~ 4 (D/50 Mpc)~^/^ /iasyr"^, which would be 
challenging to measure because of the typical lifetimes of the masers features and their 
weakness. Alternatively, a measurement of the position-velocity gradient in the low-velocity 
maser features would provide an independent estimate for their radial location within the disk 
{rsys = 0.17 [1^/50 Mpc] [A/0.27 Mpc yr'^ rad"^]-^/^ pc, where D is the distance and A = v/O 
is the velocity gradient). The limited SNR in our VLBl data precluded measurement of the 
gradient (as evident from Figure Hj). It is unclear what improvement in SNR would be 
necessary to yield a useful measurement, as there is a dearth of low- velocity features even in 
the sensitive single-dish spectra of the source. Nonetheless, the maser is time variable, and 
new spectral features may emerge with time. An independent estimate for the disk radius of 
the low-velocity features either from proper motions or position-velocity gradient could be 
used to determine a distance to NGC 3393, a result of considerable value since the galaxy is 
within the Hubble flow {vsys = 3750 kms~^) and thus might be used to es tablish a Hubbl e 



relation independent of standard candle calibrators such as Cepheids (e.g.. iGreenhilll 120041 ) . 
Over and above eventual modeling errors for VLBl data, the peculiar motion of NGC3393 (or 
the barycenter of the parent Hydra cluster) would probably limit the accuracy of inference 
for the Hubble constant from NGC3393 alone to ~ 10%. Within this, uncert ainty over the 



flow field in the vicinity of the Great Attr actor would probably dominate (IMasters et al. 



20061 ). 
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Table 1. Sources of phase error in the VLBI experiment. 



Name 


Equationt") 








(degrees) 



Uncertainty in group delay estimate 



(2ttAu) 



_3_ fTs\ 1 
27r2 \TaJ VA^i 



Residual delay error due to calibrator position error {2ttAi/) 



Residuals in atmospheric delay 

Uncertainty in fringe rate estimate 

Residual fringe rate from imperfect astrometry 

Errors in baseline length 

Errors due to imperfect calibrator astrometry 
Phase calibrator-target angular separation 



0.4 X 10"^ ■ 



B Adc 



5000 km 5 mas 



{2TTAl^)Aa^tn 



{2ttAu) t 



{2-kAu) 



0.13 X 10" 



B AOc 



5000 km 1 mas 



t 



{2-kAu) 



{2-kAu) 
K 



BAdr, 



An 



bfe^deg; 



J^Ae + ^ 

sin e 2 



0.5 - 1 

0.07 - 0.2 

< 0.6 

< 0.5 
0.3 - 1 

0.2 

0.2 

< 40 



("'Adopted from iHerrnsteinI l|l997l ) except for the last entry which is based on lCarilli &: Holdawavl l|l999l ). Au = 8 MHz is the video bandwidth (i.e., the bandwidth across 
which various calibrations are applied and over which the resulting errors are propagated), Tg ~ 1000 Jy is a representative system equivalent flux density of VLB A antennas, 
= 3 — 7Jy is the flux density of delay/rate calibr ators (i.e., 0727 -115, 0748+126, 4C 39.25), teal ~ 3.5min is the delay/rate calibrator scan duration, A6c = 0.3 — Imas is 
the uncertainty in the delay/rate calibrator position dlVIa et al.lll998h . B ~ 5000 km is an approximate baseline length, A(Tatm < 0.2nsec is the residual multiband delay after 
correcting for the zenith path delay over each antenna, Aui < 1 mHz is the residual rate after correcting for time variation in delays and fringe rates due to the troposphere and 
clock uncertainties, t = 1.5 hrs is an approximate time separation between the delay/rate calibrator scans, u>o = 22 GHz is the observing frequency, AB ~ 2 c m is a representative 
uncertainty in baseline length, c is the speed of light, A6p = 0.86 mas is the uncertainty in the phase calibrator (1055-248) position l lBeaslev et al.|[2002h . bjjni is an effective 
baseline length in km, K = 200 — 600 is a constant that depends on weather conditions (values assumed here are for the VLA site), Va < 15ms~^ is the tropospheric pattern 
speed, h ~ 2 km is the water vapor scale height, e ~ 25° is a representative antenna elevation, A9 = 2.4° is the angular separation between the phase calibrator and the target 
source, tcyc ~ 100 s is the time between successive phase calibrator observations, and Amm = 13 is the observing wavelength in mm. 

('''Phase error. On a 5000 km baseline, a phase error of 1 rad corresponds to a position error of ^ 0.08 mas. 
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Table 2. Velocities, positions, and integrated fluxes for mapped maser emission. 



Velocity^") 






Flux(^) 


(km 


1 \ 

s-i) 


(mas) 


(mas) 


(Jy) 


'^l ^4 V 


— "^1 ^7 7 


-0.3 ±0.2 


-0.9 ±0.4 


0.04 ±0.01 


^7 7- 


OlUl. 1 


-0.5 ±0.3 


-0.4 ±0.7 


0.04 ±0.02 


OlUl. 1 




-0.3 ±0.2 


-0.4 ±0.4 


0.04 ±0.01 


'^i QQ ri- 




-0.3 ±0.2 


-0.5 ±0.5 


0.04 ±0.01 




—"^90^ Q 


-0.1 ±0.2 


-1.7 ±0.4 


0.014 ±0.008 




— '^9nQ 4 


0.0 ±0.2 


-2 ± 1 


0.03 ±0.01 


'?99'? 9- 


— ^99(^ f> 


0.1 ±0.1 


-1.5 ±0.3 


0.015 ±0.006 


'^797 1- 


— ^7'^n ^ 

O 1 OU. tJ 


-0.8 ±0.2 


1.0 ±0.6 


0.012 ±0.009 


3730.5- 


-3734.0 


-1.1 ±0.3 


0.3 ±0.6 


0.02 ±0.01 


3771.6- 


-3775.0 


-1.3 ±0.2 


-0.2 ±0.3 


0.02 ±0.02 


4255.1- 


-4258.6 


— Z.O ± U.Z 


2.6 ± 0. / 


n n /I _i_ n n 1 
0.04 ± U.Oi 


4258.6- 


-4262.0 


-2.1 ±0.2 


1.3 ±0.5 


0.03 ±0.01 


4267.3- 


-4270.7 


-2.5 ±0.1 


1.9 ±0.6 


0.02 ±0.01 


4270.8- 


-4274.2 


-2.5 ±0.2 


1.5 ±0.3 


0.02 ±0.01 


4305.5- 


-4308.9 


-2.3 ±0.2 


1.2 ±0.5 


0.02 ±0.01 


4308.9- 


-4312.4 


-2.0 ±0.2 


1.0 ±0.6 


0.03 ±0.02 


4312.4- 


-4315.9 


-2.0 ±0.1 


1.2 ±0.4 


0.06 ±0.02 


4315.9- 


-4319.3 


-2.13 ±0.09 


1.6 ±0.2 


0.05 ±0.01 


4319.4- 


-4322.8 


-2.20 ±0.07 


1.6 ±0.2 


0.06 ±0.01 


4322.8- 


-4326.3 


-1.9 ±0.1 


1.8 ±0.4 


0.03 ±0.01 



Range of mapped optical heliocentric velocities. 

(''^Right ascension and declination relative to 02000 = 10'^48™'23!4660 and ^2000 = 
— 25°09'43'.'478, which lies ~ 1 mas from the estimated dynamical center (Section 4). The 
a priori position for the maser was a2oon = 10 M8'^23!45 and ^2000 = — 25°09'43'.'6 with 
uncertainty of ±0''3 from iKondratko et al.l (120061 ). 

'^'^^ Integrated flux from a fit of a two-dimensional elliptical Gaussian model to the distri- 
bution of the maser emission on the sky in each spectral channel. 
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Table 3. Fitted Line Profiles for Spectral Features Near 3880 km s~^ and 4051 km s~^ 



Xr'^"-' Velocity 


Drift 


Amplitude'^ 


FWHM(^) 


Date 


(km s-^) 


(km s"^ yr~^) 


(mJy) 


(km s~^) 






i.z oo/i.uztu.y 


Q 4- n 7 






2005 


Jan 15 






K -\- 
± Z, 


/I 4- 1 
41 It 1 


2005 


Oct 17 






1 n -I- 

±U It z 


9 4- n R 
Z.D It U.D 


2006 


Jan 22 






O -I- 

y ± z 


/I 1 4- n o 
4.i it U.y 


2006 


Mar 23 






10 4-1 

iU ± i 


O 4- n Q 

O.z ± U.o 


2006 


Apr 28 






y ± z 


6. ( ± U.y 


2006 


May 23 


oo/ 4.y ± u.z 


O.O ± u.z 


± Z 


q 1 n 7 
o.o it U. ( 


2005 


Jan 15 






41 ±2 


3.4 ±0.3 


2005 


Oct 17 






4to zn Z/ 


o.o HI U.Z/ 


2006 


Jan 22 






32 ±2 


3.2 ±0.3 


2006 


Mar 23 






29 ±2 


2.4 ±0.2 


2006 


Apr 28 






1 4- 
10 It z 


4 4- n fi 
41. Z It U.D 


2006 


May 23 


oo / o.y zt u.z 


4 4-1-00 
41.41 Zt U.Z 


1/1 K _i_ n Q 
141.0 It u.y 


K 4- n 
o.z It u.o 


2005 


Jan 15 






9"^ 4- 
ZO It o 


1.1 It u.o 


2005 


Oct 17 






1/14-9 
i4t HI z 


1 .u Hz u.o 


2006 


Jan 22 






/I 4- 9 

41: It Z 


9 4- n Q 
z.o 4z u.y 


2006 


Mar 23 






4- 1 
O It 1 


4- 1 

O Zt 1 


2006 
2006 


Apr 28 
May 23 


z.u 4uou.yi It u.u± 


u.uuo zt u.uuy 






2005 
2006 


Oct 17 
Jan 22 






119 ±5 


0.99 ±0.04 


2006 


Mar 23 






300 ±5 


0.80 ±0.01 


2006 


Apr 28 






254 ±8 


0.80 ± 0.02 


2006 


May 23 


4051.9 ±0.1 


-0.05 ±0.09 






2005 


Oct 17 






13 ±3 


1.4 ±0.4 


2006 


Jan 22 






11 ±4 


1.4 ±0.5 


2006 


Mar 23 






33 ±6 


0.7 ±0.2 


2006 


Apr 28 






23 ±4 


1.2 ±0.4 


2006 


May 23 


4053.06 ± 0.07 


-0.7 ±0.6 


32 ±5 


0.9 ±0.2 


2005 


Oct 17 
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Table 3 — Continued 



Xlt^") Velocity(^) Drift('=) 


Amplitude^'^) FWHM('^) 


Date 


(km s~^) (km yr~^) 


(mJy) 


(km s ^) 






7±3 


1.4 ±0.8 


2006 Jan 22 








2006 Mar 23 




6±2 


1.4 ±0.7 


2006 Apr 28 








2006 Ma.v 23 



Reduced chi-squared for a least squares solution that determines velocity drifts, ampli- 
tudes, peak centroids, and widths of all Gaussian components at all epochs simultaneously. 

(''^Velocity of a Gaussian component at the reference epoch (2005 January 15 and 2005 
October 17 for ~ 3880 km s~^ and ~ 4051 km s~^ features, respectively). 

('^^Velocity drift of a Gaussian component. 

Amplitudes of a Gaussian component at the listed epochs. 

(^^FuU widths at half maximum of Gaussian components at the listed epochs. 



0.10 
0.08 
0.06 — 
0.04 
0.02 



0.00 






3000 



3500 4000 
Helio Opticol Velocity (km s"') 



4500 







3100 3150 3200 

Hello Opticol Velocity (km s"') 



3750 3800 3850 
Helio Opticol Velocity (km s"' 



4250 4300 4350 

Helio Opticol Velocity (km s"') 



Fig. 1. — Top: A weighted average of five total-power spectra of the maser in NGC3393 
obtained with the GBT on 2005 October 17, 2006 January 22, 2006 March 23, 2006 April 28, 
and 2006 May 23 with an effective resolution of 108kHz (~ 1.5kms~^) and la = 0.6 mJy. 
Bottom: Comparison of VLBI imaged power and total power spectrum. The total-power 
spectrum was obtained with the GBT on 2005 January 15, with an effective resolution of 
108 kHz (~ 1.5kms~^). Angle-integrated flux density (filled circles) and la measurement 
uncertainty (shaded boxes) are shown for each 3.4kms~^ spectral channel with detectable 
emission. Though separated by ~ 9 months, the single-dish and VLBI flux measurements 
agree in general to ~ la, the exception being the apparent absence of 3880 km s^^ 
(marked by the arrow and for which we report a velocity drift of 5 ± lkms~^ yr" 
time of the VLBI observations. 



emission 
^) at the 
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I I I I I 1 1 1 1 

1 0-1-2 -3 -4 

Relative RA (mas) 



Fig. 2. — Distribution of maser emission in the nuclear region of NGC3393. Position uncer- 
tainties are la, and the colors of the maser spots indicate heliocentric optical line-of-sight 
velocity in accordance with the bar on the right. The dotted line in the color bar shows 
the adopted systemic velocity of 3750 km s~^. The adopted location for the dynamical cen- 
ter (black circle) is the weighted mean for the low-velocity maser features. A line fitted 
to the distribution of maser emission on the sky (P.A.~ —34°) is close to orthogonal to 
the kpc-scale radio jet (black arrows: P.A.~ 45°, Morganti et al. 1999; ~ 56°, Schmitt 
et al. 2001b ) and to the axis of the NLR (dashed cone : P.A.~ 55° with an opening an- 
gle of ~ 90° ISchmitt fc Kinney! Il996l : ICooke et al.ll2000l ). The coordinates are relative to 
^2000 = 10M8'"23!4660 and (Jaooo = -25°09'43'.'478. At a distance of 50Mpc, 0.24 pc sub- 
tends 1 mas. 



- 21 - 




3870 3875 3880 3885 3890 4046 4048 4050 4052 4054 4056 4058 4060 

Optical Heliocentric Velocity (km s"') Optical Heliocentric Velocity (km s"') 



Fig. 3. — Maser spectra of features at ~ 3880 km s~^ and ~ 4051 km s~^ obtained with 
the Green Bank Telescope and the results of Gaussian component decomposition. The 
feature near 4051 km s~^ was not detected on 15 January 2005, emission in the immediate 
vicinity of the systemic velocity, ~ 3750 km s~^, was too weak to allow Gaussian component 
decomposition and is not shown. 
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4500 - 




-4 



2 2 

Impact Parameter (mas) 
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Fig. 4. — Position- velocity diagram for the detected maser emission in NGC 3393. A fit of 
the Keplerian rotation law to the high-velocity features (filled grey curves) yields a mass of 
(3.1 ± 0.2) X 10^ Mq enclosed within 0.36 ± 0.02 pc (1.48 ± 0.06 mas). From the measured 
velocity drift (a = 5 ± Ikms^^yr^^), we infer the radius of the systemic feature at ~ 

^GMBu/a = 0.17 ± 0.02pc, 

0.005 ± 0.001 radyr-i 



3880 km s ^ (not detected at the VLBI epoch itself )of Vgys 
which corresponds to Qsys = 



J /^sys 



a^ly{GMBH] 



1/4 



'GMbh, 

(nearly vertical, filled grey cone). The systemic velocity estimates are cited by the NED (3(T 
uncertainties are plotted). In our calc ulations, we adopted t he systemic velocity from 21 cm 
line measurements of 3750 ± 5kms~^ ( ITheureau et al.lll998h . 



